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1. Gauge-fixing as a sampling procedure 

Gauge theories have a very interesting structure. Take a theory of a set of fields 0. If a fixed 
field value at each space-time point is given, this defines a configuration O [0] of the fields. In a 
gauge theory there exist local transformations — ^ + 5 of the fields (/) such that for the set 
of configurations called an orbit, reachable by these transformations the corresponding path 
integral remains invariant, i. e. all correlation functions have the same value. Especially, this implies 
that all correlation functions not invariant under a gauge transformation vanish. This vanishing is 
the realization of Elitzur's theorem [1]. If this is the only symmetry of theory, the non- vanishing 
correlation functions define the non-trivial set of observables of the theory. 

In principle, their calculation is all that is necessary to determine all experimental conse- 
quences of a theory. One possibility to do this are, e. g., lattice calculations. However, not all 
interesting cases can be solved with such methods efficiently. This has lead to the development of 
methods using a different approach, and which include, e. g., perturbation theory and functional 
methods. 

In these cases, one removes the gauge symmetry, i. e., breaks it explicitly, in such a way as 
that any obsei-vable is not modified. This is done by instead of just integrating over the whole orbit 
in the path integral with a flat weight, integrating over it with some suitable non-flat weight: A 
non-trivial sampling of the orbit is introduced. This includes the extreme case of a 5-f unction as 
weight to pick out a single element of the orbit, a single gauge copy. In general, if the sampling 
includes only a subset of the orbit, this subset will be called here the residual gauge orbit. If the 
residual orbit contains gauge copies which cannot be deformed into each other by infinitesimal 
gauge transformations, these are called Gribov copies [2]. 

For the implementation of this gauge-fixing procedure several possibilities exist. In lattice 
calculations it is indeed possible to perform this sampling explicitly [3]. In perturbation theory, 
this can usually be achieved by the inclusion of auxiliary fields, the so-called ghost-fields [4]. The 
situation becomes complicated for non-lattice non-perturbative methods, like functional methods 
[3, 5]. After the discussion of examples of such weighting procedures in section ^, the current state 
of the art in the continuum will be revised in section ^ An interesting possibility is to include 
residual symmetries in the weighting process, which will be outlined in section ^. Finally, a few 
concluding marks will highlight the current challenges in section ||. For the sake of being explicit, 
only Yang-Mills theory will be discussed here, though most of this can be found easily in other 
theories as well, e. g. QCD. 



2. Residual symmetries 

A starting point to illustrate the concept of residual symmetries are the covariant gauges in 



2 



Local and global gauge-fixing 



Axel Maas 



perturbation theory [6], described by the (Euclidean) path integral 



j m^&c&cff{A^,c,c)e-^'''''^^ (2.1) 



Df = 5'"d,+gf^'A'^^, 

where the gauge fields A and the ghost and anti-ghost fields c and c interact with a coupling strength 
g and live in the adjoint representation of the gauge algebra codified by the structure constants 
jrabc^ and is a parameter which modifies the sampling of the orbit with the Gaussian weight 
factor. The ghost field contribution only compensates, due to their origin from a Jacobian [6], the 
sampling in a way as to keep the values of the observables invariant. Thus, the last two terms 



of the Lagrangian ( |2.2[ ) represents the perturbative sampling procedure, the gauge-fixing. As a 
consequence, con^elation functions not invariant under local gauge transformation no longer vanish. 

These gauges have several continuous residual global symmetries, under which the sampling 
procedure is still invariant. One are global color rotations. As a consequence, any correlation 
function having a color direction still vanishes. To change this would require a further reduction 
of the sampled orbit. This is usually done when adding scalar fields in 't Hooft-type gauges [6], 
though this is not necessary [7]. Similarly, there is a trivial scale symmetry of the ghost fields, 
associated with ghost number conservation [6]. 

Far more interesting is the third symmetry, the BRST symmetry. This global symmetry trans- 
forms both ghost and gluon fields [6]. Since it modifies the gluon fields without altering any gauge- 
invariant observables, this is necessarily a gauge transformation, though one which also alters the 
ghost fields. It thus connects the gauge copies, and the change in the ghost fields is necessary to 
transform the weight accordingly. 

In the case of the perturbative Landau gauge, i. e. the limit — )■ 0, the BRST symmetry can no 
longer act non-trivially on the gluon fields, as there is one and only one gauge copy per orbit which 
satisfies the perturbative Landau gauge condition d^A"^ = 0. However, there are still non-trivial 
ghost transformations, which leads to an additional ghost- anti-ghost symmetry [8], essentially ro- 
tating ghost and anti-ghost into each other. 

The situation becomes vastly more interesting when going non-perturbative. Then Gribov 
copies [2, 9] appear. Thus BRST transformation can now act once more non-trivial even in Landau 
gauge, and one regains a non-trivial and non-perturbative BRST symmetry [7]. However, this 
requires, as it likely seems [10-12], to include all Gribov copies with a flat weight. This leads to 
a number of complications, which can only be appreciated after investigating the appearance of 
Gribov copies a bit more closely. 



3. Landau-gauge examples on the lattice 

Once non-perturbative, there are many Gribov copies satisfying the Landau gauge. The so 
obtained residual gauge orbit turns out to have a highly non-trivial structure [3, 13]. In particular, 
it is possible to define a bounded and convex Gribov region D., enclosed by the so-called Gribov 
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horizon, in which all eigenvalues of the Faddeev-Popov operator —d^^Lf^ we strictly positive and 
only one vanishing on the horizon. Since every orbit passes through this region [14], it is possible 
to restrict the sampling of the orbit further to this region by the introduction of a ©-function on the 
Faddeev-Popov operator, with 0(0) = 1. 

Still, many Gribov copies remain inside this Gribov region. To deal with them, several differ- 
ent gauge conditions have been invented [3]. These either attempt to identify a single Gribov copy 
as the representative of the gauge orbit, being the absolute Landau gauge [15] and its inverse [16], 
based on the so-called fundamental modular region of minimal gauge field norm or by externaliz- 
ing either correlation functions [17] or eigenvalues [18] of the Faddeev-Popov operator. Since such 
constructions are notoriously complicated to construct in the continuum, this will not be discussed 
further here, see [3] for more details. 

It appears more interesting to pursue gauges averaging over the Gribov region with some 
prescription [3, 5]. In general, any such averaging procedure is performed in the same way, i. e.. 



by rewriting the perturbative expression (2.1 ) as [3, 19] 



<^>= lim / ^Au^c^c€?(A„,c,c)0(-5uD"*)e--^^'^-^«w(A„,c,c), 

where w is an appropriately chosen weight functions, which includes a normalization such that any 
observable remains unchanged. There are a few prominent examples in use. 

In lattice calculations, a convenient choice is w = 1, i. e. averaging over the first Gribov region 
with a flat weight. Since usual Landau-gauge-fixing algorithms appear to identify Gribov copies 
with equal probability, this is in practical calculations simplified to take a single, random represen- 
tative for each configuration, leading to the minimal Landau gauge [3]. In continuum calculation 
another choice is based on the (yet unproven) assumption that there exists a weight function such 
that @{—dfxD"^)w{Afji,c,c) = 5(— 5^D^^) [13], i. e. averaging with a flat weight over the horizon 
only. This choice has the advantage that it can be rewritten as a local Lagrangian using further 
auxiliary ghost fields. 

In contrast to these ai^e gauges which average with some weight, either exponential [3, 5] 
or Gaussian [19], with some ai^gument over the Gribov region. Of course, once appropriately 
normalized, this does not change observables. Such an averaging includes a control parameter 
A, essentially the width of the sampling function. This is a second gauge parameter, besides the 
perturbative one t, . 

To be concrete, consider the following two possibilities [3, 19] 

wi = exp + ^ I d''xd''yd;c''{x)dy{y)^ (3.1) 

W2 = exp - ^ I d'^xA^Al^ , (3.2) 

where the are appropriately chosen normalizations. The first case [3, 5] is a boundary term, 
essentially driving the ghost dressing function at zero momentum to a desired value with Ai being 
a Lagrange parameter. Especially, the value A, = is a fixed point returning the minimal Landau 
gauge, while the values — oo and +oo drive the ghost dressing function at zero momentum to the 
minimum and maximum possible values. In the second case [19], the Gribov copies are weighted 
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with the norm of the gauge field. When A2 — > °°, this will put all weight on the Gribov copies in 
the fundamental modular region, i. e. the region of least gauge field norm. Again, the fixed-point 
value A2 = returns the minimal Landau gauge. 

To visualize the impact of both prescriptions, the gluon propagator, the ghost dressing function, 
and the running coupling in three dimensions are shown for several values of the As in figure [I] [20]. 
First of all, it is clearly visible how, within statistical errors, the cases for A, = coincide with the 
minimal Landau gauge, as expected from the general arguments. 



Concerning the weight function (p.l[), there are several observations, in agreement with earlier 
results [5]. When changing the gauge pai'ameter A, away from zero towards large positive and 
negative values, the results move towards those of the extremal maxB and minB gauges developed 
in [17], respectively. In these single-copy gauges the ghost propagator was extremalized. This is in 
agreement with the results and reasoning of [5]. In contrast to original hopes [17], but in agreement 
with insights based on the BRST construction in [3, 21, 22], no qualitative change is observed in 
the infrared region. Especially the running coupling, though within large errors, remains vanishing. 
However, this will require more detailed investigations in the future for a more definite statement 
[20], though any change would be surprising. 



The situation for the weight function ( |3.2[ ) is similar. Indeed, as argued [19], the results tend 
to agree better with the absolute Landau gauge with an increasing value of the gauge parameter A2, 
though the dependence is rather weak. Nonetheless, also this gauge fulfills in its lattice implemen- 
tation its intended purpose. 

However, in both cases it should be kept in mind that not all Gribov copies are guaranteed to be 
found in lattice simulations [5]. Thus the observe results are something of a lower limit to what the 
impact of Gribov copies really is. Thus an additional systematic uncertainty remains. Especially, it 



is never possible using the weighting function (3.2) in a lattice calculation using any contemporary 



algorithm to achieve a better approximation of the absolute Landau gauge than by any other method 
[3, 16, 23-25] of searching the absolute minimum, since it can never be better than the best Gribov 
copy found. But this does not exclude the possibility that new algorithms based on this idea maybe 
developed which can improve this. 

4. Taking the continuum limit 

Besides the technical problem of finding all Gribov copies, the problem remains of how to 
incorporate weights like ( [3.1| ) and (3^) into continuum calculations. The 0-function is known 



not to change the form of functional equations [26]. This implies [5] that these equations contain 
the solutions both inside, outside, and in the whole of gauge field configuration space. These 
appear to consequently yield more solutions than observed on the lattice in the first Gribov region 
[21, 27]. It remains thus an unsolved problem how to classify the continuum solutions of functional 
methods according to whether they belong to the first Gribov region or not, otherwise than by the 
(necessarily approximate) comparison to lattice results. Even the approach of [13] cannot guarantee 
this, as it has not been ruled out that a similar replacement is justified from outside the first Gribov 
region. 



Concerning the further weight (3.1), the weight function is actually a surface term. Such 



terms turn into boundary conditions [3, 28] in functional equations, which indeed appear to yield 
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Gluon propagator with type 1 weighting 



Gluon propagator with type 2 weighting 
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Figure 1: Preliminary results on the gluon propagator (top panels), ghost dressing function (middle panel), 



and running coupling (bottom panel) for the weight functions (3.1) (left panels) and (3.2) (right panels) 
for various values of the gauge parameters A,-, compared to both the minimal Landau gauge and the absolute 
Landau gauge. Results are from three dimensions with V = (7.9 fm)^ (42^ lattice) at a = 0.189 fm (j8 = 3.92) 
[20]. Note the caveats of [5] on the problem of finding Gribov copies in lattice calculations. More details 
will be available elsewhere [20]. 
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precisely the same results as the use of the weight function in lattice calculations [3, 21]. However, 
a formal understanding, or even proof, is still lacking. 

The second weight (3.2) has so far only been used in the context of a perturbative calculations 
[19], similar- in spirit to the approach of [13]. Since the term is not a boundary term, it is not 
equivalent to a genuine boundary condition, though some arguments exist [17, 23], how this may 
be achieved indirectly. 

In total, it remains still an open challenge how to formally correctly implement the same gauge 
condition on the lattice and in the continuum, though the rather good quantitative agreement ob- 
tained between both approaches [3, 5, 21], as well as to other results, is already quite encouraging. 



5. Concluding remarks 

In summary, it becomes clear that the structure of the gauge orbit and how to sample it, is 
irrelevant for any observable. This is, as it must be: The introduction of the gauge symmetry is 
only a technical tool to obtain a local field theory, which otherwise can only be described in temis 
of non-local objects, e. g. Wilson lines on a finite lattice. This also implies that any valid way 
of sampling the gauge orbit is equally admissible, and it is valid as long as it leaves observables 
invariant. 

While this is a well-developed subject in perturbation theory [4], it only becomes now more 
and more relevant beyond perturbation theory. The reason is simple. Various non-perturbative 
methods have reached a degree of maturity which makes them very capable of investigating physics. 
However, many of them, especially in the continuum, use a gauge-fixed framework, just like per- 
turbation theory. Furthermore, all non-perturbative methods employ some kind of approximation, 
which cannot be systematically controlled. E. g. on a lattice this maybe the continuum limit, which 
is usually only addressed by extrapolations, and in continuum methods it may be the truncation of 
equation hierarchies [3]. Thus, it is desirable to compare the results of different methods at the most 
elementary steps, to improve systematic reliability. Since these are often now gauge-fixed quanti- 
ties, it is necessary to guarantee for a comparison that gauge orbits are sampled in the same way, 
since otherwise there maybe fundamental differences between the results. This has led historically 
to a number of problems [3, 8]. 

Besides the problems induced in this way by the sampling procedure, there is also an advan- 
tage. Choosing a suitable sampling may make problems tractable or simpler [3, 7]. This has been 
widely used in perturbation theory [6]. Hai^nessing this possibility in non-perturbative calcula- 
tions, e. g. by formulations like ( ^ ) and (3.2), will possibly help to solve many non-perturbative 
problems more efficiently. 
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